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2016  AFOSR  Final  Report 

Victor  S.  Batista,  Clifford  P.  Kubiak  and  Tianquan  Lian 

Research  progress. 

We,  the  collaborative  team  of  Batista  (Yale  University),  Kubiak  (University  of 
California,  San  Diego)  and  Lian  (Emory  University),  have  made  significant  progress  in 
the  understanding  of  the  operation  of  catalysts  for  CO2  reduction  on  surfaces,  including 
Re(I)  bipyridyl  complexes,  Ni(cyclam)  complexes,  and  pyridine-based  molecules  as  well 
as  development  of  in  situ  spectroscopic  techniques  for  studying  catalytic  reactions  during 
the  last  grant  period.  Each  project  involves  an  intimate  combination  of  synthetic, 
electrochemical,  computational,  and  spectroscopic  techniques,  especially  as  related  to 
sum  frequency  generation  spectroscopy  (SFG)  and  in  situ  electrochemical  SFG. 


1.  Coupled  experimental  and  theoretical  determination  of  orientation  of  Re(I) 
bipyridyl  complexes  for  CO2  reduction  on  Au  and  on  TiCh. 

We  have  recently  determined  the  binding  orientations  of  Re(I)  bipyridyl 
complexes  with  different  anchoring  groups  on  two  different  surfaces.  First,  we  used  sum 
frequency  generation  spectroscopy  and  calculations  of  SFG  spectra,  based  on  density 
functional  theory,  to  investigate  Re(R-2,2'-bipyridine)(CO)3Cl  (R  =  4-cyano  or  4,4'- 
dicyano)  electrocatalysts  binding  to  gold  electrodes.  Understanding  the  orientation  of 
these  complexes  is  critical  for  redox  state  transitions  and  catalytic  turnover. 


Figure  1.  Left.  DFT  optimized  monodentate  geometry  for  the  dicyano  Re  complex  with 
its  axial  CO  facing  the  surface  with  a  tilt  angle  of  63°  relative  to  the  surface  normal. 
Right.  Representative  PPP-polarized  SFG  spectra  of  monolayers  of  the  dicyano  Re 
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complex  adsorbed  onto  gold  thin  films  (black  circles)  with  DFT  simulations  of  the  SFG 
spectra  using  the  geometry  on  the  left  (red). 

We  found  that  the  electrocatalysts  lean  on  the  Au  surface  and  orient  the  plane  of 

the  bipyridine  ligand  at  63°  relative  to  the  surface  normal  by  computing  several  binding 
geometries  in  DFT  and  determine  which  gives  orientations  and  thus  SFG  spectra  that  best 
match  the  experimentally  measured  spectra.  These  findings,  recently  published  in  the 
Journal  of  Physical  Chemistry  C  (Clark  2016),  demonstrate  the  capabilities  of  the 
approach  including  rigorous  spectroscopic  and  theoretical  methods  for  revealing  the 
conformation  and  orientation  of  CO2  reduction  catalysts  bound  to  electrode  surfaces.  This 
information  has  consequences  for  catalytic  turnover,  as  the  labile  chloride  ligand  and 
therefore  reactive  site  of  the  metal  is  oriented  away  from  the  surface,  thus  no  crowding 
will  occur.  Current  efforts  are  focused  on  exploring  alternative  anchoring  groups 
including  thiol,  phosphonate,  and  isocyanide  as  well  as  other  substrate  surfaces  including 
Si. 

Our  other  effort,  recently  published  in  the  Journal  of  Physical  Chemistry  C  (Ge 
2016),  was  focused  on  surface-induced  anisotropic  binding  of  the  same  catalyst,  but  with 
carboxyl  anchoring  groups  (ReCOA)  on  rutile  single-crystalline  Ti02  (110)  surfaces 
relative  to  (001)  surfaces  using  vibrational  sum  frequency  generation  (SFG) 
spectroscopy.  The  SFG  signal  shows  an  isotropic  distribution  with  the  rotation  of  the 
TiC>2  (001)  surface  relative  to  the  incident  plane,  but  an  anisotropic  distribution  with  the 
rotation  of  the  TiC>2(110)  surface  (Figure  2,  top).  By  combining  these  results  with  ab 
initio  SFG  simulations  and  with  modeling  of  ReCOA-TiCk  cluster  binding  structures  at 
the  density  functional  theory  level,  we  revealed  that  the  origin  of  the  optical  anisotropy 
for  ReCOA  on  the  TiCk  (110)  surface  is  associated  with  the  binding  preference  of  ReCOA 
along  the  [-110]  axis.  Along  this  direction,  the  binding  structure  is  energetically  favorable 
due  to  the  formation  of  proper  hydrogen  bonding  between  the  carboxylate  group  and 
passivating  water  molecules  on  the  Ti02  (110)  surface  (Figure  2,  bottom).  The  tilt  angle, 
defined  by  the  bipyridyl  ring  angle  relative  to  the  surface  normal,  of  the  catalyst  is  found 
to  be  a  combination  of  26°  and  18°,  which  corresponds  to  an  aggregate  at  high  surface 
coverage  with  minimal  steric  interactions.  Since  molecular  orientation  in  these  systems 
can  significantly  influence  their  catalytic  performance,  the  ability  to  control  the  molecular 
ordering  through  careful  selection  of  the  surface  structure  and  symmetry  would  be 
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extremely  beneficial.  This  investigation  of  catalytic  molecules  on  surfaces  with  different 
symmetries  helps  to  elucidate  how  the  degree  of  molecular  ordering  depends  on  surface 
structure  and  therefore  is  of  great  practical  interest,  which  will  inform  our  selection  of 
substrates  in  future  work. 
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Figure  2.  Top.  Polar  plots  of  the  azimuthal  dependence  of  the  amplitude  for  (a)  the  a'(l) 
mode  of  ReCOA  on  TiCk  (001)  and  (b)  a'(l)  and  a'(2)  modes  of  ReCOA  on  TiCk  (110). 
Black  squares  (a'(l))  and  red  circles  (a'(2))  are  experimental  results;  Solid  lines  are 
theoretical  results.  The  amplitude  of  the  a'(l)  mode  at  <j)  =  90°  for  ReCOA/TiCk  (110) 
system  is  normalized  to  one  with  all  other  data  points  for  each  mode  scaled  accordingly. 
Bottom.  Binding  structures  of  two  ReCOA  complexes  as  a  dimer  on  the  TiC>2  (110) 
surface  along  the  [-110]  axis  from  the  side  (left)  and  from  above,  with  the  ReCOA 
complexes  represented  by  green  rectangles  indicating  the  binding  sites  (right).  Hydrogen 
bonding  to  the  ReCOA  complexes  is  depicted  using  dashed  bonds.  The  atoms  are  colored 
as  follows:  with  green  =  Cl,  silver  =  Ti,  cerulean  =  Re,  white  =  H,  red  =  O,  gray  =  C,  and 
blue  =  N.  The  tilt  angles  of  each  catalyst  are  indicated  by  the  angle  formed  by  the  bpy 
ring  (black  dotted  line)  and  the  surface  normal  (orange  dotted  line).  The  black  arrows 
indicate  the  crystal  directions  on  the  surface  for  each  panel. 
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2.  The  surface  as  a  ligand:  the  electronic  origins  of  catalytic  enhancement  of 
Ni(cyclam)  activity  by  amalgated  gold  electrodes  (Au/Hg). 

Ni(cyclam)  has  been  the  focus  of  many  investigations  as  a  CO2  reduction  catalyst, 
but  it  suffers  from  poisoning  from  CO.  However,  a  dramatic  enhancement  of 
electrochemical  C02  conversion  to  CO,  catalyzed  by  [Ni(cyclam)](PF6)2,  is  observed  on 
Au/Hg  electrodes  when  compared  to  other  metallic  electrodes  such  as  Zn.  Thus,  our 
recent  efforts,  in  an  article  being  prepared  for  submission,  have  focused  on  understanding 
the  underlying  reaction  mechanism  responsible  for  the  increase  in  turnover  frequency 
using  DFT  calculations  on  cluster  models  with  or  without  CO  bound.  The  computed 
trans- III  conformer  structures  on  the  Au/Hg  and  Zn  electrodes  are  given  in  Figure  3 
where  similar  effects  are  found  between  conformers.  We  found  that  that  Hg  provides 
favorable  interactions  for  enhanced  reaction  kinetics  by  stabilizing  the  CO  unbound  form 
of  the  catalyst,  predominantly  with  the  trans- III  conformation,  due  to  reduced  Ni-CO 
reverse  dative  interactions  when  compared  to  the  complex  in  solution,  or  bound  to  other 
metallic  surfaces.  The  higher  CO  stretching  frequency,  shorter  CO  bond,  longer  Ni-C(O) 
bond,  more  linear  Ni-C-0  angle,  and  spin-density  closer  to  the  formal  oxidation  number 
of  Ni(I)  correlate  with  lower  CO  desorption  free  energies.  These  findings  are  particularly 
relevant  to  the  design  of  electrode  surfaces  for  activation  of  electrocatalytic  transition 
metal  complexes  with  common  surface  interactions  that  could  improve  catalyst 
performance  under  cell  operating  conditions. 
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Figure  3.  DFT  optimized  configurations  of  [Ni(cyclam)]+  on  Au/Hg  (top)  and  Zn 
(bottom),  without  CO  (left)  and  CO  bound  (right).  Color  key:  yellow:  Au,  silver:  Hg, 
blue:  N,  white:  H,  green:  Ni,  red:  O,  dark  gray:  carbon,  purple:  Zn. 

3.  Proton-coupled  electron  transfer  (PCET)  in  the  reduction  of  4,4’-bipyridine  with 
CdSe  quantum  dots  and  in  the  reduction  of  pyridinium  on  platinum  for  the 
conversion  of  CO2  to  methanol. 

Our  efforts  into  the  understanding  of  PCET  in  pyridine-based  molecules  are  based 
on  two  fronts:  reduction  of  bipyridine  by  CdSe  quantum  dots  and  of  pyridinium  by  Pt 
electrodes.  In  the  first  case,  recently  published  in  the  Journal  of  the  American  Chemical 
Society  (Chen  2016),  we  studied  the  photo-reduction  of  4,4’-bipyridium  (bPYD)  using 
CdSe  quantum  dots  (QDs)  as  a  model  system  for  interfacial  PCET.  We  observed  ultrafast 
photo-induced  PCET  from  CdSe  QDs  to  form  doubly  protonated  [bPYDH2]+"  radical 
cations  at  low  pH  (4-6).  Through  studies  of  the  dependence  of  PCET  rates  on  isotope 
substitution,  pH  and  bPYD  concentration,  the  radical  formation  mechanism  was 
identified  to  be  a  sequential  interfacial  electron  and  proton  transfer  (ET/PT)  process 
(Figure  4a)  with  a  rate  limiting  pH  independent  electron  transfer  rate  constant,  kmi,  of 
1.05  ±  0.13  x  10 10  s'1  between  a  QD  and  an  adsorbed  singly  protonated  [bPYDH]+.  In  the 
presence  of  sacrificial  electron  donors,  this  system  was  shown  to  be  capable  of  generating 
[bPYDH2]+"  radical  cations  under  continuous  illumination  at  405  nm  with  a  steady-state 
photoreduction  quantum  yield  of  1.1  ±  0.1%  at  pH  4.  DFT  geometry  optimizations  of 
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4,4 ’-bipyridine  on  a  model  CdSe  cluster  were  carried  out  in  the  presence  of  bound  3- 
mercaptopropionic  acid  (MPA)  and  were  compared  to  equivalent  calculations  of 
methylviologen.  The  binding  geometry  is  given  in  Figures  4b  and  4c.  The  optimized 
structure  for  bPYD  shows  asymmetrical  hydrogen  bonding  so  that  only  one  H+  is  bonded 
to  the  bPYD  (1.22  A)  while  the  other  H+ is  merely  hydrogen  bonding  to  it  (1.59  A), 
corroborating  the  role  of  [bPYDH]+  in  electron  injection.  The  hydrogen  bonding 
interaction  with  the  MPA  prevents  bPYD  from  getting  close  to  the  surface,  as  opposed  to 
methylviologen  with  its  permanent  positive  charge,  suggesting  that  the  faster  injection 
into  methylviologen  is  due  to  its  ability  to  bind  closer  to  the  surface.  Theoretical  studies 
of  the  adsorption  of  [bPYDH]+  and  methylviologen  on  QD  surfaces  revealed  important 
effects  of  hydrogen  bonding  with  the  capping  ligand  (3-mercaptopropionic  acid)  on 
binding  geometry  and  interfacial  PCET.  The  mechanism  of  bPYD  photo-reduction 
reported  in  this  work  may  provide  useful  insights  into  the  catalytic  roles  of  pyridine  and 
pyridine  derivatives  in  the  electrochemical  and  photoelectrochemical  reduction  of  CO2. 
This  work  established  that  QDs  can  be  used  as  a  model  system  for  studying  interfacial 
PCET  mechanisms. 
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Figure  4.  Interfacial  proton-coupled  electron  transfer  from  QDs.  a)  Scheme  of  proton- 
coupled  electron  transfer  from  QDs  to  bipyridine,  b)  and  c)  Structure  of  bPYD  on  the 
CdSe  {100}  cluster  model  from  above  (panel  b)  and  from  the  side  (panel  c).  The  color 
code  is  as  follows:  white  =  H,  red  =  O,  gray  =  C,  yellow  =  S,  blue  =  N,  brown  =  Se,  and 
tan  =  Cd. 
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Figure  5.  Top:  proposed  mechanism  of  electrochemical  reduction  of  CO2  on  Pt  (1 1 1)  in 
an  aqueous  imidazole  solution  (pH  =  5.2).  Bottom:  Thermodynamic  cycle  (ii.l)-(ii.3) 
used  to  obtain  the  free  energy  change  due  to  reduction  of  imidazolium  adsorbed  to  the  Pt 
surface  to  form  Pt-hydride  on  the  surface. 

In  the  second  case,  as  published  recently  in  collaboration  with  Bocarsly  and  co¬ 
workers  in  Topics  in  Catalysis  (Liao  2015),  we  have  computationally  found  that  the 
mechanism  of  CO2  reduction  by  imidazolium  is  similar  to  the  reduction  by  pyridinum,  as 
given  in  Figure  5,  top.  It  was  found  to  form  H  adsorbates  on  the  electrode  surface  through 
a  one-electron  PCET  reaction  (-0.68  V  vs.  SCE  as  computed  by  Figure  5,  bottom)  rather 
than  a  one-electron  reduction  to  the  imidazolium  radical.  The  cycle  includes  surface 
desorption  of  imidazolium  to  dihydrogen  and  imidazole  in  the  aqueous  solution  (ii.2),  and 
dissociative  desorption  of  H  to  the  Pt  surface  (ii.3).  The  reduction  potential  of 
imidazolium  to  H  and  imidazole  (ii.2)  was  predicted  to  be  -0.60  V,  in  good  agreement 
with  the  reduction  potential  observed  by  CV  experiments,  while  the  free  energy  changes 
for  reactions  (ii.l)  and  (ii.3)  approximately  cancel  each  other.  These  results  are 
particularly  valuable  since  they  suggest  that  H  adsorbates,  as  recently  corroborated 
through  our  other  recent  collaborative  work  with  isotope  effects  on  electrochemistry, 
published  in  the  Journal  of  the  Electrochemistry  Society  (Zeitler  2015),  could  be 
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generated  at  low  overpotentials  by  either  of  the  two  electrocatalysts.  The  H  adsorbates 
might  react  as  Pt-hydrides  and  reduce  CO2  through  a  proton  coupled  hydride  transfer 
(PCHT)  mechanism,  as  recently  proposed  for  the  electrocatalytic  reduction  of  CO2  to 
formic  acid  in  the  presence  of  pyridinium  where  the  electrophilic  attack  of  CO2  onto  the 
surface  hydride  is  activated  by  the  Bronsted  acid  in  solution. 

4.  Vibrational  relaxation  dynamics  of  Re(R-2,2'-bipyridine)(CO)3Cl  (R  =  4,4'- 
dicyano)  electrocatalyst. 

We  have  started  to  investigate  how  metal  electrodes  and  interfacial  solvent 
environment  affect  vibrational  energy  relaxation  dynamics  of  catalysts,  a  process  that  is 
intimately  coupled  to  chemical  reactions.  We  have  studied  vibrational  relaxation 
dynamics  of  Re(R-2,2'-bipyridine)(CO)3Cl  (R  =  4,4'-dicyano)  electrocatalyst  on  gold  and 
Ti02  surfaces  and  in  different  solvents.  Vibrational  dynamics  of  Re  complex  monolayer 
bound  to  the  gold  surface  was  probed  by  IR  pump/SFG  probe  spectroscopy  (Figure  6a). 
Vibrational  relaxation  dynamics  of  Re  complexes  adsorbed  on  nanoporous  TiCk  thin  film 
and  in  solution  were  measured  by  IR  pump/IR  probe  transient  absorption  spectroscopy. 
Our  initial  results  showed  that  the  both  metal  electrodes  and  solvent  facilitates  the 
vibrational  relaxation  of  CO  stretching  modes  in  the  Re  complex  (Figure  6b).  Ongoing 
theoretical  calculations  are  aimed  at  revealing  the  mechanisms  by  which  the  substrate  and 
solvent  affect  the  vibrational  relaxation  of  Re(I)  electrocatalysts  using  anharmonic 
coupling  constants  as  yielded  by  electronic  structure  calculations. 
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Figure  6.  a)  Normalized  time  resolved-SFG  difference  spectra  of  three  CO  stretching 
modes  of  Re(I)  electrocatalyst  on  Au  at  indicated  averaged  delay  times.  Solid  circles 
represent  experimental  data;  solid  lines  represent  fits,  b)  Kinetics  of  the  GSB  of  in-phase 
symmetric  stretch  mode  for  Re(I)  electrocatalyst  on  Au,  Ti02,  and  in  different  solvents. 
Symbols  represent  experimental  data;  lines  represent  fits. 

5.  In  situ  electrochemical  SFG:  electric  field  effects  at  the  electrode. 

Our  latest  effort  involves  investigating  the  electrochemical  Stark  effect  via  SFG 
spectroscopy  in  order  to  correlate  dependence  of  vibrational  frequencies  of  chemisorbates 
on  the  electrode  potential  to  useful  information  about  the  double-layer  structure  at  the 
electrode/solution  interface.  We  are  investigating  the  Stark  effect  from  the 
nitrile/isocyanide-terminated  self-assembled  monolayers  (SAMs)  on  the  gold  electrode 
surface.  A  series  of  molecules  with  different  terminal  groups,  symmetry,  or  length  was 
chosen  to  systematically  study  the  local  electrical  fields  near  the  electrode  surface.  To 
further  understand  the  Stark  effect  on  the  SAMs  systems,  theoretical  calculations  were 
performed  to  simulate  the  vibrational  frequency  shift  of  nitrile  or  isocyanide  groups  in 
response  to  local  electrical  fields.  For  example,  Fig.  7a  shows  a  scheme  of  1,4-phenylene 
diisocyanide  (PDI)  adsorbed  on  the  electrode  surface.  Fig.  7b  shows  the  potential- 
dependent  SFG  signals  of  isocyanide  groups.  At  0  V  vs.  Ag/AgCl,  two  peaks  around 
2120  (free  -NC)  and  2170  cm-1  (Au  bonded  -NC)  were  observed.  At  more  negative 
potentials,  a  blue  shift  of  the  free  N=C  mode  and  a  red  shift  of  the  Au  bonded  N=C  mode 
were  observed  (Fig.  7c).  Such  a  difference  could  be  explained  by  the  different  responses 
to  the  electric  fields  for  these  two  groups,  which  have  opposite  dipoles  (Fig.  7d). 
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Figure  7.  (a)  Scheme  of  PDI  adsorbed  on  gold  electrode  surface,  (b)  Potential-dependent 
SFG  spectra  of  the  PDI  SAM.  Circles  represent  experimental  data,  solid  lines  represent 
fitting  results,  (c)  Experimental  peak  positions  of  free  -NC  and  Au  bonded  -NC  under 
different  potentials,  (d)  Calculated  peak  positions  of  free  -NC  and  Au  bonded  -NC  under 
different  electric  fields. 
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